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Abstract 

It was previously found that in a certain parameter subspace of scalar- 
tensor theories emerging from massive gravity, the only stable field configura- 
tion created by static spherically symmetric sources was one with cosmological 
asymptotics. Moreover, these backgrounds were shown to be sub-luminal ev- 
erywhere in the space; in contrast to the common believe that these theories 
are necessarily super luminal in the vicinity of a static source. In this work 
we complete that analysis by extending it to cover the whole parameter space 
of these scalar-tensor theories. We find that the stability argument renders 
the asymptotically flat backgrounds unrealizable, forcing once again for cos- 
mological asymptotics. In the case of pressureless sources these backgrounds 
are stable. However, they get destabilized in the presence of positive pressure, 
larger than a critical density. Even on the self-accelerated background, on 
which the scalar mode decouples from sources, in the region occupied by the 
source it acquires an elliptic equation of motion. Therefore, we conclude that 
the only parameter space which is not ruled out, by solar system measure- 
ments, is the one considered in Berezhiani et al. (arXiv: 1302.0549), namely 
the one for which the scalar and tensor modes can be diagonalized via local 
transformations. 

We also reinvestigate the scale at which perturbation theory breaks down 
in a general Galileon theory. We show that the Vainshtein mechanism success- 
fully redresses the strong coupling scale to a small one, just like in the cubic 
Galileon, despite the cancellations occurring in the special spherically sym- 
metric case. We emphasize that even if these tests were performed at scales 
at which perturbation theory broke down, these could not be interpreted as a 
lower bound for the graviton mass. 



1 Introduction and Summary 



There exists a two-parameter family of theories of massive gravity that propagate 
five degrees of freedom in four dimensions [1, 2]. In the decoupling limit, it gives rise 
to a fascinating class of scalar-tensor theories. There is a one-parameter sub-class 
of these theories, for which the scalar mode can be completely decoupled from the 
tensor using the invertible field redefinition h^ u — > h^ u + ni]^ + ad^Tid v Ti / A3, [1]. As 
a result of this de-mixing, the longitudinal mode of the graviton becomes described 
by the so-called quartic 'Galileon' model [3] 1 , supplemented with a novel coupling 
to the matter stress-tensor d fJ ,7rd u 7rT IIU , where 7r denotes the helicity-0 mode of the 
graviton and T^ u denotes the matter energy-momentum tensor, [1]. After partial 
diagonalization, the Lagrangian describing the scalar sector reduces to 

where A 3 = (Mpim 2 ) 1 / 3 is the strong coupling scale and a is the free real parameter 
of this sub-class of theories. 

The thorough study of the spherically symmetric configurations [7] (for related 
work on spherical symmetric solutions see [9]) showed that the stability of the spher- 
ically symmetric configurations in the presence of the static source forces the free 
parameter a to be positive. Otherwise, the last term of (1) gives rise to a ghost-like 
kinetic term in the high density regions. Furthermore, it was shown in [7] that, for 
a > 0, the theory does not admit asymptotically flat classical solutions. Instead, 
the screened it configuration at short distances (within the Vainshtein region [8], see 
also [10] for a recent review on the Vainhstein mechanism) matches a cosmological 
background at large distances. We emphasize that this matching effect is not related 
to the last term of (1) and would be present even in its absence. In addition, it was 
explicitly shown that these backgrounds are subluminal, as opposed to the common 
believe that in massive gravity the configurations recovering GR at short distances 
necessarily exhibit superluminal propagation 2 . 

In the present work we would like to complete the analysis of [7] by extending 
it to the full parameter space of the decoupling limit of massive gravity. In general, 
there is a nonlinear mixing term between the helicity-0 and -2 modes which cannot 
be removed by a local field redefinition. This makes the assay slightly more involved 
and gives rise to a qualitative change in the conclusions. Namely, for generic pa- 
rameters of the theory, we show that asymptotically flat backgrounds created by 



1 This type of theories were first discovered in [4] in the context of the DGP scenario [5], see 
also [6] 

2 We emphasize however that the 'issue' of superluminality in massive gravity has not (yet) 
been connected to that of acausality in any rigorous way. Configurations on which closed time-like 
curves could form seem to live beyond the regime of validity of the effective theory, [11]. 
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static spherically symmetric sources exhibit a gradient instability. This should be 
contrasted with the case considered in [7]; where asymptotically flat configurations 
were unstable as well, although the instability was ghost-like. 

The presence of this gradient instability in the general case forces us to give 
up asymptotic flatness, and we instead focus on the self-accelerated background 
configuration. The latter was already shown to be interesting for its unique property 
to decouple the n mode from matter (at leading order), [13, 14], not to mention 
the obvious phenomenological significance of such backgrounds. The absence of a 
direct conformal coupling tcT on top of this self-accelerating solution implies that no 
Vainshtein mechanism or other screening effect ought to be invoked. Furthermore, 
after using Einstein's equations, this decoupling leads to the following contributions 
to the kinetic term of the scalar mode in the vicinity of the source 



The examination of (2) shows that, in case of a dust-like source, both contributions 
to the kinetic term 7r 2 are healthy; hence, the scalar mode is not a ghost. However, 
it is straightforward to see that in case of the positive pressure the second term 
of (2) gives rise to a negative contribution to the gradient energy. Moreover, if 
the pressure is larger than the critical density p c = M P iA|, something which is 
common for astrophysical sources or indeed even for the atmosphere! this negative 
contribution overcomes the positive one, making the background unstable. This 
argument can be extended to other asymptotic cosmologies, as long as we recover 
General Relativity at short distances the same instability will arise. 

Therefore, if the scalar-tensor theory considered here were to reproduce all the 
important effects of the full theory, then we would conclude that the only phe- 
nomenologically viable theory of massive gravity is the one with (3 = and a > 0. 
However, the scalar-tensor theory may not be necessarily capturing all the impor- 
tant properties of massive gravity. In particular, the presence of a nonzero time-like 
component of the helicity-1 field (A ^ 0) may give rise to a new class of spherically 
symmetric solutions not captured by the scalar-tensor sector. 

In the second part of this manuscript, we investigate the strong coupling issue 
in theories such as massive gravity and Galileons. The fact that the Vainshtein 
mechanism relies on irrelevant operators (in the Wilsonian sense) to be important 
raises the question of the control one has over perturbation theory. As previously 
investigated, [4] over a background configuration for the scalar field n = n + 5n, 
the strong coupling scale gets redressed, symbolically A re d re ssed ~ A3(<9 2 7ro/A3) a , 
with a power a depending on the exact model. In the cubic Galileon that arises 
in DGP, a = 1/2, while in models with higher order Galileon interactions, one 
could have a = 2/3. This redressing allows to raise the strong coupling scale to 
~ (1cm) -1 in the cubic Galileon case, however this scale is still low enough to 
wonder what happens at energy scales above that. Furthermore in the quartic 
Galileon and some models of massive gravity, it has been argued that the redressing 
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only raises the strong coupling scale to ~ (0.4km) _1 making the theory fully non- 
perturbative at extremely low energy scales, [16]. We show here that the difference 
between the cubic Galileon and the quartic Galileon/massive gravity is actually not 
as pronounced as previously found, and the redressed strong coupling scale in the 
quartic Galileon/massive gravity is actually closer to ~ (30cm)" 1 , i.e. at most one 
order of magnitude below the cubic Galileon. 

Putting this subtlety aside, one can raise the interesting question of what happens 
when reaching energy scales higher than the redressed strong coupling scale, and 
perturbation theory runs out of control. Following Vainshtein's original argument 
(rather than its specific realization within a Galileon theory), we do expect the effect 
of the graviton mass and in particular the effect of the helicity-0 mode to become 
smaller and smaller as one gets to higher and higher energies. We emphasize however 
that there is no sense in which one could use the break of perturbation theory as a 
lower bound for the graviton mass as has been done in the literature [16]. 

The paper is organized as follows. Section 2 is dedicated to the description of 
the framework. In section 3 we analyze the spherically symmetric configuration 
and study its stability both for asymptotically flat backgrounds and for the self- 
accelerated one. We also discuss the stability for more general cosmologies. Finally, 
we reexamine the strong coupling issue in section 4, and open the discussion for a 
deeper understanding in section 5. 

As we were finalizing this paper, it came to our attention that the related work 
was being conducted in [12]; which has some overlap with our work. 

2 The Theory 

The scalar-tensor sector of massive gravity, in the decoupling limit, is described by 
the following Lagrangian density 




where, we have denoted the helicity-±2,0 modes by h^ v and n respectively. 

The three identically conserved symmetric tensors X^fll] depend on second 
derivatives of the helicity-0 field Tl^ = d^d u Ti in the following way, 

x (1) = --e apa e 13 n r 

v( 2 ) — I £ a n £ TT JT 

V"( 3 ) — £ a PT £ ^TT aTT TT c 

The Lagrangian (3) is invariant under the linear diffeomorphisms Sh^ = + ctC/i 
up to a total derivative, while being exactly invariant under the global Galilean 

symmetry 5n = v^x^. 
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Recently, it has been shown [15] that the decoupling limit of massive gravity ex- 
hibits the same non-renormalization properties as the Galileon theory [3, 4]. Namely, 
the coefficients a and j3 do not get radiatively corrected, within the effective theory. 

In [7], it was found that the f3 = parameter subspace of (3) does not possess 
a stable asymptotically flat solutions in the presence of the spherically symmetric 
static sources. Instead, the stable screened ir configuration inside the Vainshtein 
region is smoothly matched to the cosmological backgrounds (with various equations 
of state) at large distances. In this work, we would like to perform the similar analysis 
for the rest of the parameter space (i.e. j3 ^ 0), which is qualitatively different from 
its [3 — counterpart. In particular, in the present case the mixing between the 
different helicity states cannot be undone by means of a local field redefinition. This 
structure is so far characteristic to massive gravity [2] and has not been observed in 
other modifications of General Relativity. 



3 Static Spherically- Symmetric Configurations 

Throughout this work we consider a star-like source of finite size R and uniform 
density p. Without loss of generality, the static spherically symmetric configuration 
can be found by assuming the following ansatz for the metric perturbations around 
a Minkowski space-time in spherical coordinates 

hoo = a(r); hij = f(r)6ij. (4) 

The most general ansatz (up to Galilean transformations) for the helicity-0 mode, 
which leads to the static spherically symmetric metric (4) (see appendix A), is given 
by 

n = ^Alt 2 + n (r). (5) 

After integrating Einstein's equations, i.e. the equations for h^, once and using 
vanishing initial conditions in the origin we obtain 

rf' = -^ + Ay\(l-a\-2^), (6) 

ra' = -^~ r + A^r 2 (c - A(l + 2ac) - 6f3c\ 2 - 2/3A 3 ) . (7) 

Following the same procedure for the longitudinal mode we arrive at the following 
equation 

3(1 + 2ac)X - Q(a + a 2 c - 4/3c)A 2 + 2(a 2 - 4/3 - 20af3c)\ 3 - 60/3 2 cA 4 - 12/3 2 A 5 
2 (^) 3 (1 + 2ac + 12/3cA + 6/3A 2 ) + c Outside the source 
2 (if) 3 ( X + 2ac + 12 / 3cA + 6 ( 3X2 ) + c Inside the source 
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where we have defined 



x < , ( M \ 1/3 

being the Vainshtein radius and ir' = <9 r 7r (r). The classical backgrounds in the 
time independent case (c = 0) have been previously studied in [17, 18]. It has been 
established that for generic parameters (8) possesses two types of solutions inside the 
Vainshtein region. In particular, since the factor (r*/r) 3 is large at short distances, 
(8) requires either (i) A > 1 or (ii) the vanishing coefficient of (r*/r) 3 . Moreover, 
it has been shown that the case (i) corresponds to a screened gravitational field at 
short distances, hence contradicting empirical data. 

In this work we study the general class of spherically symmetric configurations 
(5), created by static sources, against instabilities. 

3.1 Asymptotic Flatness 

As it has been already mentioned for generic a and j3 the equation of motion for 
the longitudinal mode has two types of solutions at short distances, r<Cr t . One of 
these solutions is 

a = -/r 1/3 -, (10) 

r 

which corresponds to a screened gravitational field, i. e. a field for which the New- 
tonian 1/r behavior is absent, at short distances even for nontrivial c; as it is easy 
to see from (6) and (7). We dismiss this solution for obvious reasons 3 . 

The other solution corresponds to the case when the coefficient of (r*/r) 3 in (8) 
vanishes on the background, that is 

1 + 2ac + 12/3cA + 6/3 A 2 = 0. (11) 

Which evidently gives us a solution, with A = const; hence, recovering GR at short 
distances with high precision (see (6) and (7); the normal 1/r behavior of gravity 
at large distances precisely cancels). This serves as a motivation for studying the 
stability of the generic background with constant A. Within the Vainshtein region 
the leading contribution to the gradient energy, to the second order in fluctuations 
reads as follows 

£(2) = -3/3( c + A) (^) 3 [2(d r 5n) 2 - (d n 5n) 2 ] . (12) 

Notice the relative minus sign between the radial and angular terms 4 . As a result, 
we deduce that the only way to avoid the gradient instability is to have A ~ — c at 

3 Moreover, for c = this solution matches the cosmological background, as it was noted in [17]. 
4 Here, we would like to emphasize that (12) comes from the WXffi term of (3); the other 
contributions to the quadratic gradient term are subdominant. 
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short distances. This condition, combined with (11) leads to 

_ a ± y/g 2 + 6/3 



(13) 



If we further require asymptotic flatness, we are led to the following system of 
conditions 

A 00 (l-aA 00 -2/3A 2 ) = 0, (14) 

c - Aoo (1 + 2ac) - 6/?cA^ - 2/3A^ = 0, (15) 

where A^ denotes the value of A at spatial infinity. It is easy to show that these 
conditions possess a nontrivial solution only in the parameter space 

a 2 

fi = ~Y- (16) 

However, according to [20], in this parameter subspace the flat vacuum with A = — c 
is infinitely strongly coupled. 

Hence, the analysis of this section leads us to the conclusion that the theory 
under consideration does not possess a stable asymptotically flat solution sourced 
by a static and spherically symmetric source. 



3.2 Asymptotically de Sitter 

Having failed to introduce a source on flat space, it is natural to ask what happens 
when we substitute the flat space by a self-induced de Sitter vacuum. We look for 
the de Sitter space solution in static slicing for which the ansatz adopted in the 
previous section is applicable here as well. 

Therefore, we are looking for the spherically symmetric solutions to (8) with de 
Sitter asymptotic. This means that at large distances the effective energy momentum 
tensor, coming from the mass term of the graviton, will have the equation of state 
of the cosmological constant. As a result, using (6)-(8), we arrive at 

-a± ^/a 2 + 6(3 
c =-A(r^oo) = v — . (17) 

Then, the eq.(8) at arbitrary distances factorizes and takes the following form 

(A + c)P(A) = 12/3(^) 3 (A + c) 2 , (18) 

where, P(X) is the polynomial of the fourth order; which does not vanish when 
A = — c, unless /3 — 0, /3 — —a 2 / 6 or (3 = — a 2 /8. The first possibility leads to 
a ghost-like instability of the de Sitter space [13], in the second case the helicity-0 
loses its kinetic term, while the third one necessarily leads to an asymptotically 
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flat background (which has already been discarded due to infinite strong coupling). 
This implies that the only solution for the generic parameters (except the above- 
mentioned ones) is the one with the trivial (source independent) n profile around 
the source 

n = ~Alx^. (19) 

This may be traced to the fact that, according to [13], on de Sitter space the helicity-0 
mode does not have the kinetic mixing with the helicity-±2. Namely, the Lagrangian 
density re-calculated on the self-accelerated background is given by 

£ = ~\^StK, + «(a, P)^nD, + 

-^^^(dnynn + + ^T, v . (20) 

Here, and it denote the deviations from the de Sitter background for the tensor 
and scalar modes respectively. The curvature of the background is given by H 2 oc m 2 
and k is some function of the parameters. The exact expressions for them are 
inessential for current discussion and in case of interest the reader is directed to 
[13]. 

Again, the only asymptotically decaying solution to the equations of motion, 
following from (20), is ir = 0. While h^ u satisfies the linearized Einstein's field 
equations 

G - = w, T >- (21 > 

The problem in case of certain sources arises from the third term of (20), which by 
integration by parts can be rewritten in the following form 

^P-G^d^. (22) 

In case of the sources smaller than their Vainshtein radius this term will overwhelm 
the kinetic term inside the sources, precisely like in [7] . Hence the 7r becomes a ghost 
unless k > 0; surprisingly, this is exactly what we need in order to have the healthy 
kinetic term outside of the source as well (20). 

Therefore, we conclude that the theory admits the presence of the arbitrarily 
dense objects on de Sitter space. However, this is so only for pressure-less sources; 
otherwise an extra caution is in order. Namely, if the source has a pressure, larger 
than the critical density, then the dominant contribution to the quadratic term comes 
entirely from (22). Moreover, in case of the positive pressure this contribution has 
a wrong sign (as it follows from (21)) in front of the gradient term, leading to the 
gradient instability of the configuration. Unfortunately, most of the localized sources 
have sufficiently large pressure to realize this instability; even the atmosphere around 
the earth has the pressure 10 14 times larger than the critical density. 
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3.3 Other Cosmologies 

We would like to start this section by pointing out that, up to this point the situation 
is qualitatively similar to the /3 = case, see [7]. Namely, the presence of natural 
sources seems to be destabilizing the asymptotically flat, as well as asymptotically 
self-accelerated spaces (in case of f3 — there was no stable self-accelerated space 
to begin with, see [13]). However, in that case the situation has changed once we 
allowed for cosmological backgrounds with nontrivial equation of state [7]. In this 
section we would like to explore this possibility for the /3 ^ parameter space. 

First, we reiterate that the stability of the background forces upon us the con- 
dition A ~ — c within the Vainshtein region (see eq.(12)). As a result, the leading 
contribution to the quadratic (in perturbations) Lagrangian vanishes. In order to 
find the next order, no n- vanishing result, we have to find the corrections to the 
background itself. The correction to A is parameterized as A = — c + 5X; where, as it 
is easy to see from (8), 5 A is of order of (r/r*) 3 . While the expression for the metric 
degrees of freedom at next to leading order is obtained by substituting A = — c to 
the right hand side of (6) and (7); leading to the correction ~ A 3 r 2 . In other words, 
the background next to the leading order has the following form 



A ~ -c + O 



(23) 



A/f A 3 r 2 

a ~ - A 3 r 2 c (-1 -ac + 2/3c 2 ). (25) 

By simple analysis of these corrections we conclude that the leading contribution to 
kinetic terms of fluctuations comes from the background 

A = -c, (26) 

while for / and a we should take into account the corrections present in (24) and 
(25). To this end, the background relevant for the quadratic terms is similar to the 
one for de Sitter space (see the previous subsection). With the distinction that in 
the case of de Sitter space the background (26) was exact, while in the current case it 
is merely the leading expression for the background. Despite this minor distinction, 
since (26) is the relevant piece of the background at short distances, the conclusions 
of the previous section apply here as well. Namely, the background will exhibit the 
instability either inside the source or right outside of it. It should also be mentioned 
that the mixing between the scalar and tensor modes of the graviton are irrelevant, 
since to the leading order we will have the pure n kinetic term, while the hir mixing 
appears only in the sub-leading approximation. 

We would like to conclude by noting that in case of j3 <C a there seems to be 
another branch of solutions at short distances. Namely, one could think of having f3 
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so small that the right-hand side of (8) could have been neglected at short distances; 
such a small values for the parameter are technically natural since it does not get 
renormalized by quantum loops [15], so this would not represent a fine-tuning of 
the parameters. For simplicity, let us concentrate on c = and a = 1 as it is quite 
straightforward to generalize the argument. This branch of solutions requires A ^> 1, 
which under the assumption that the left hand side is negligible gives 

2A 3 ~ 12/3 2 A 5 A ~ ±-^. (27) 

V6/3 

However, in order this assumption about the right hand side to be valid, we need to 
make sure that 

1 /r*\3 

"TvT > 



(s) ■ (28 > 



.R, 

Taking into account (6), (7) and the expressions given above, it is easy to see that 
this branch of solutions gives an unacceptably large deviation from GR. Hence, it is 
ruled out on phenomenological grounds. 

4 Strong Coupling for general Galileons 

In the previous section, we have shown that for massive gravity, stability of the 
spherically symmetric solutions imposes a non-flat, i.e. cosmological, asymptotic 
behaviour. This is however not necessarily the case for a more general Galileon 
theory where the coefficients between the different operators are relaxed and there 
is no mixing with the helicity-2 mode of the graviton. In this section we investigate 
the redressing of the strong coupling scale on spherically symmetric configurations 
and discuss the small departure from spherical symmetry. We start with a generic 
Galileon setup and compute the redressing of the strong coupling scale both in the 
cubic and the quartic Galileon. Within the Vainshtein regime, the stable solution in 
massive gravity covered by a > is essentially a special case of the quartic Galileon 
(although the asymptotics outside the Vainshtein region are different). 

4.1 Asymptotically Flat Spherically Symmetric Galileons 

To start with, we consider a quartic Galileon Theory in four dimensions 

^/^£xiM + ^ r ' <29) 

where T is the trace of the stress-energy tensor of any external sources and the 
Galileon Lagrangians are given by 

C 2 = (<9tt) 2 (30) 
£ 3 = {dTifUu (31) 
A = (dn) 2 ((Dtt) 2 - (d M cU) 2 ) . (32) 
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In massive gravity (with /3 — 0), — 1, C3 = —a and C4 = a 2 /3, but we leave 
the C3 and C4 coefficients arbitrary for now and simply set — 1. We are mainly 
interested in the Vainshtein screening of this Galileon due to the Earth, so as first 
approximation, the background source may be considered as spherically symmetric, 
T = — M e S^(r) + 5T, leading to a background configuration for the Galileon field 
which is also spherically symmetric 



n(t,r) = 7T (r) + — 0(t,f) , 
where the background field satisfies the simple algebraic equation, 

2 o / / / w 3 



+ 



2C3 



+ 



2c 4 
AT 



*0 \ r , 



12tt Mpi r 3 



(33) 



(34) 



If the quartic Galileon is present, C4 7^ 0, then close to the source, (at small r), the 
last term on the left hand side dominates and the background solution takes the 
form 

1/3 



A 2 



1 M a 



+ 0(A 3 r) 



while if c 4 



127rc 4 M P1/ 

0, the background field acquires a different profile, 



A 2 



M ff 



1/2 



+ 0(A 3 r) 



(35) 



(36) 



127rc 3 M P1 A 3 r x 

The perturbations around this background configuration evolve with the following 
kinetic matrix Z^, 



= --Z^{r)d,<pd u <P + 



with 



12 



= 1 + 4c3?r o + 6C4^ 



Ai r 



V) = 



LA 

3r 2 dr 
1 d 



A^ r' 

1 + ^ + 
A3 r 



/2 



2r dr 



18c 4 7To- 
4C3 7Tq 6C4 7Tq 2 



Af r 



A3 1 r 2 



(37) 

(38) 
(39) 
(40) 



When the quartic Galileon is present, one can see immediately that the leading con- 
tribution to the angular direction vanishes, Z u ~ Z rr ~ (A3r)~ 2 (M e /Mpi) 2 ^ 3 while 
Znn ~ (A 3 r) _1 (M e /Mpi) 1//2 . This leads to a few subtleties in the quartic Galileon 
but as we will see later, the resulting redressed strong coupling scale is nevertheless 
barely affected by this subtlety. The reason for that is that the same cancelation 
responsible for the hierarchy Zqq <C Z tt is also responsible for cancelling the leading 
contribution to the operator that would naively arise at the lowest energy scale. As 
a result the redressed strong coupling scale is larger than naively anticipated. 



11 



4.2 Cubic Galileon 



4.2.1 Redressing from the Earth 

We start by focusing on the cubic Galileon and set c 3 = 1/3 for simplicity. In that 
case 

z tt ~ z„. ~ z nn ~ (^|-) ^-L^ ee Z ffi . (41) 

The canonically normalized field is then 

(42) 

leading to the cubic interaction 

4 3) = xs (^) 2 n ^ = (^) 2 n< ^ = xr (^) 2 n ^ • ( 4 3) 

So the redressed coupling scale due to the screening of the Earth is given by 

/ M \ 1/4 1 

A e = A 3V / ^~A 3 — — —H7- (44) 

® 3V e 3 V47rMpJ (A 3 r) 3/4 1 ' 



(A 3 rf 

Taking the strong coupling scale to be that associated with an infrared modified 
theory such as DGP or massive gravity for which A 3 = (ifgMpi) 1 / 3 where H ~ 
10~ 33 eV, we have A 3 ~ (1120km) -1 , and so at the surface of the Earth the redressed 
scale is (as previously found in [19]) 

A ffi ~ 10 7 A 3 ~ (4 cm)" 1 . (45) 
4.2.2 Redressing from Nearby Sources 

When testing the Newton's law using torsion balance at submillimeter scales, the 
experiment itself and nearby sources can further contribute to the screening of the 
Galileon field. For concreteness, let us consider a source of mass M e (which could 
symbolize the experiment itself or a nearby source) localised a distance p e from 
the core of the experiment. The coupling to that new source leads to a new field 
configuration vr e (p) on the top of that 7r (r), given by 

where p is the distance from the source M e . Going through the same analysis 
presented previously, the redressing of the strong coupling scale due to the mass M e 
is then 

/ M \ 1/4 1 

A -Mw^J (47) 
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As a possible example, if we consider that the local effects could be mimicked by a 
mass of 10kg localized 1cm away from the center of the experiment then, 

A e ^ 4 A e . (48) 

In itself this is not a huge enhancement, but it simply serves to show that as one 
goes within the apparatus itself, its different components help raising the scale, and 
screening the force. 



4.3 Quartic Galileon - Massive Gravity 

In massive gravity, when (3 = 0, the decoupling limit resembles that of a quartic 
Galileon with an additional coupling to matter. As we will discuss in section 5, in 
that case, the only relevant scale in the decoupling limit is not A 3 = (m 2 Mpi) 1//3 but 
rather 

Since the coefficient a does not renormalize [15], a can in principle depart signif- 
icantly from unity, hence disentangling the strong coupling scale A which appears 
in the decoupling limit, with the graviton mass m. In what follows we will thus set 
c 4 = a 2 /3 as is the case in massive gravity and work with the scale A which can in 
principle be independent from the graviton mass. 



4.3.1 Angular Subtleties 

As mentioned previously, in the purely spherically symmetric case, when in the 
vacuum Zqq <C Z rr ~ Z tt , which leads to a few subtleties in treating this case. In 
what follows we consider the case where 

6 ( 1 lid \ 

and Z nn = e Z ffi , (51) 

with e < 1. 

• If we only consider the effect from the Earth itself, and are interested in the 
redressed scale in the vacuum outside the Earth, then as seen earlier, 

1 AR * fcwxW-»f 7 =A™V (52) 



\ Mpl J 

• However, the Earth itself is not perfectly spherically symmetric, and just taking 
into account the flatness of the Earth, (which is of the order of 5 ~ 0.0033), we would 
have instead more realistically e ~ 5 2 ~ 10 -5 . Furthermore, the presence of other 
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sources near the experiment itself will completely break the symmetry and more 
realistically, we would expect e ~ 1 at the level of the experiment itself. However 
for consistency we keep e as an arbitrary parameter for now, with e < 1. 

Since the gradient along the angular direction are not redressed with the same 
scale as along the radial direction, we first need to rescale the angular directions as, 
[16] 

(t,r,e, V ) = (if,e 1 / 2 6,e 1 / 2 0), (53) 

this is not something one could do globally, but if we are only interested in what 
happens in a small region of space near say an experiment, and not for all angles, 
the rescaling can be done locally. The kinetic term is then of the form 

J d 4 x£ = J d 4 x(-±Z'»(r)d li <i>d v <i> + --?) ~ J d A x^~dJd^+--^j , (54) 

with = y/eZ^ (j). In terms of the canonically normalized field 0, the quartic 
and cubic Galileon leads to cubic interactions of the form (focusing only on the 
interactions that arise at the lowest energy scale), 



/ 



d 4 *£ mt d / d* x [ ^mMdU) + + ^ir ) (W + ■■■) (55) 




v A e,i A e,2 



with the redressed interaction scales, 



1/3 



Ae,i = A 2 ^J (^e) V2 (57) 

Ae,2 = A 2 ^y /3 (eZe) 1/2 ~A6 1/3 ( e Z ffi ) 1/2 . (58) 

Now going back into the original coordinates, if we read the redressed scale as a 
scale in the orthoradial direction then, 

A©,i,2 = -4A e)li2 (59) 
v e 

This leads to the following redressed interaction scales 

A ffi ,i ~ AZ'J 3 (60) 

A e , 2 = Ae 1 / 3 ^. (61) 

Being extremely conservative and taking the naive value for e as given in (52), we 
get 

A e> i ~ 5 x 10 6 A ~ (20cm)" 1 (62) 

A e , 2 ~ 3 x 10 5 A ~ (30cm)" 1 , (63) 
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assuming A 3 ~ (H^Mpi). This differs from the result of (0.4km) 1 obtained in [16] 
by about three orders of magnitude. The reason for the discrepancy lies in the fact 
that the operator (<9^0) 2 is not enhanced by the coefficients as is the case for 
the other operators of that form. This is due to the specific structure of the Galileon 
interactions. As we show here, this operator has a coefficient going as vr^r) rather 

n' (r) 

than - ai - i . Since the hierarchy between Zqq and Z tt was precisely coming from 

the fact that around spherically symmetric configuration ^'(r) <C ^p- and so that 
operator comes in with a much larger energy scale than one could have anticipated 
at first sight 5 . 

Being even more conservative and translating instead the redressed scales as 
distance scales along the radial direction, we would then have instead A^i^ ~ A-®,i,2, 
leading to a smaller scale (or larger distance scale along the radial direction). Then 
for a more realistic value of e < 10 -5 , we would get instead A ffiil ~ (70m) _1 , for 
A ~ (HqMpi) 1 / 3 as a radial distance scale (and A ffi2 ~ (lm) _1 ). However any 
realistic kind of matter located around the experiment itself would completely wash 
out this hierarchy and set e close to 1, giving back Ae^^ ~ few x (10cm)" 1 for 
A ~ (HqMpi) 1 / 3 . Since any test of Newton's law requires an apparatus which itself 
is of the order of a fraction of a metre, and is itself relatively massive, we expect the 
mass present at these distance scales to redress the strong coupling scale further. 
More work needs to be done to fully explore the influence of other environmental 
factors to the redressing of this scale. 

5 The Bigger Picture 

As we have seen in the previous sections, massive gravity and generic Galileons have 
a relatively low energy scale at which perturbation theory breaks down. However, 
it is not because perturbation theory breaks down that we suddenly expect large 
corrections to Newton's law. If anything, we would expect the Vainshtein mechanism 
to work even better in this non-perturbative regime and to decouple the field even 
more. The scales A 3 , A e or A re d reS sed m &y or may not be interpreted as a cutoff in a 
usual sense, i.e. it is not necessarily the case that new degrees of freedom come in 
at this scale. One may imagine it is the scale of strong coupling, analogous to Aq CD , 
and non-perturbative tools must be developed to understand what happens at these 

5 When dealing with the mixing h^ u X^ v , a confusion of similar nature was made in [16]. Indeed, 
thanks to the very precise nature of the interactions in massive gravity, operators of the form 
(9 2 7r) 3 always appear in such a combination so that they form a total derivative, or in other words, 
there are no operators of the form (<9 2 7r) 3 in massive gravity, which is the essence of its ghost- 
free construction [1, 2]. Instead when conside ring sm opcrj itor of the form h(d 2 ir) 3 , its relevant 
contribution arises after integrations by parts to give rise to the operator (d 2 h)(d7r) 2 (d 2 Tr), which is 
of course similar to *i? A " yQ/3 9 M 7r9 a 7r-D 1/ D^7r, this is an operator of dimension 7 rather than 9 which 
becomes relevant at a larger energy scale. However more importantly, as mentioned previously, 
there are no stable solutions which exhibit the Vainshtein mechanism when this mixing term is 
present. 
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scales. For this reason it is important to stress that the scale at which perturbation 
theory breaks down cannot be used to put a lower bound on the graviton mass as 
has been argued in [16]. 

Interestingly for massive gravity the stability analysis of this paper constrains 
us to the restricted Galileons, considered in [7], for which /3 — 0. In this case the 
only scale entering into the decoupling limit is A = A^/a 1 ^ 3 . Since both a and (3 are 
free parameters which satisfy a non-renormalization theorem [15], the scale of strong 
coupling is completely independent of the mass of the graviton. Furthermore since 
essentially all observational constraints so far, constrain only the scale A = As/a 1 / 3 
6 and not A3 directly, they do not impose any constraints directly on the graviton 
mass. The later must be constrained by cosmology where the decoupling limit is 
not appropriate, rather than solar system/astrophysical gravitational tests [21, 22]. 

The results of this and the previous paper [7] have important implications for the 
discussion on the possible existence of superluminalities, and UV completion. For 
instance, as argued in [23], generic Galileons violate the conditions for analyticity of 
the S-matrix in Minkowski space-time. However, the very definition of an S-matrix 
assumes the switching off of interactions at infinity consistent with the assumptions 
of asymptotic flatness. What we have seen is that introducing a single source into 
the theory is in conflict with asymptotically flat boundary conditions, at least in 
the decoupling limit, for stability reasons. As such the Minkowski S-matrix for 
the decoupling limit theory is not an appropriate description, and the existence 
or not of its analyticity is a mute point. This is not to say that these theories 
do not need to have a fundamentally different non-perturbative/UV completion 
than traditional Wilsonian effective field theories, but rather that their failure to 
satisfy usual analyticity properties in the decoupling limit by no means precludes 
the existence of a non-perturbative completion, even one for which the physics is 
potentially fundamentally (sub-)luminal. 
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Appendix A 



In order to find the most general ansatz for 7r, one has to satisfy the following 
conditions: 

(i) Spherical symmetry requires ir(r,t) to be a function of r and t only. 

(ii) For the configuration to be static the effective energy- momentum tensor 
(T^5 = G^ v ) must satisfy the following 

T f = and d t T^ = 0. (A-I) 

Here, the first expression requires the vanishing momentum density while the latter 
one assures time-independence of the energy density. One could notice that because 
of the energy-momentum conservation (Bianchi's identity), these two constraints 
are not independent from each other. However, nothing prohibits requiring both of 
them individually. Notice, that these conditions are the direct consequence of (4). 

It follows from the Einstein's equation, that the condition of vanishing momentum- 
density reads as follows 

-i + ^n fcfc + M(n» fc -n„n„)) % 



a 6/3 . 



n 0j = o. (a-ii) 



For simplicity let us imagine that U 0i and are independent variables. Since, in 
that case we have a linear system of algebraic equations at hand, with U Q j unknown. 
There exists a non-trivial solution to (A-II) if and only if the determinant of the ma- 
trix in brackets vanishes. Otherwise, we end up with dod r 7i = 0, which is consistent 
with our final ansatz (A-V). With the assumption of spherical symmetry 



7r' , „ 7T 



U ij = — 6 H + niUjin" - — ), (A-III) 



the above-mentioned condition becomes 




£ - =0- (A-IV) 



From (A-IV) and the time independence of follows that the most general ansatz 
relevant to us is 

7r = T(t) + 7r (r), VT,vr . (A-V) 

If we further require the time-independence of the effective stress tensor T^ ff , we 
arrive at the following ansatz 

7r = |A 3 t 2 + 7r (r), (A- VI) 

where we could add terms constant and linear in time, however those terms are 
irrelevant because of the Galilean symmetry. 
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